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Abstract 

Casimir energy is calculated for the 5D electromagnetism and 5D scalar theory 
in the warped geometry. It is compared with the flat case (arXiv: 080 1 . 3064) . A new 
regularization, called sphere lattice regularization^ is taken. In the integration over 
the 5D space, we introduce two boundary curves (IR-surface and UV-surface) based 
on the minimal area principle. It is a direct realization of the geometrical approach to 
the renormalization group. The regularized configuration is closed-string like. We do 
not take the KK-expansion approach. Instead, the position/momentum propagator 
is exploited, combined with the heat-kernel method. All expressions are closed-form 
(not KK-expanded form) . The generated P/M propagators are introduced. Rigorous 
quantities are only treated (non-perturbative treatment). The properly regularized 
form of Casimir energy, is expressed in a closed form. We numerically evaluate A(4D 
UV-cutoff), a;(5D bulk curvature, warp parameter) and T(extra space IR parameter) 
dependence of the Casimir energy. We present two new ideas in order to define the 
5D QFT: 1) the summation (integral) region over the 5D space is restricted by two 
minimal surfaces (IR-surface, UV-surface) ; or 2) we introduce a weight function 
and require the dominant contribution, in the summation, is given by the minimal 
surface. Based on these, 5D Casimir energy is finitely obtained after the proper 
renormalization procedure. The warp parameter oo suffers from the renormalization 
effect. The IR parameter T does not. In relation to characterizing the dominant 
path, we classify all paths (minimal surface curves) in AdSs space. We examine the 
meaning of the weight function and finally reach a new definition of the Casimir energy 
where the 4D momenta ( or coordinates) are quantized with the extra coordinate as 
the Euclidean time (inverse temperature). We comment on the cosmological constant 
term and present an answer to the problem at the end. Dirac's large number naturally 
appears. 

PACS: PACS NO: 04.50.+h, ll.lO.Kk, 11.25.Mj, 12.10.-g 11.30.Er, 
Keywords: Warped geometry, position/momentum propagator, heat-kernel, Casimir 
energy, sphere lattice, Renormalization of boundary parameters, Cosmological constant 
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Figure 1: The configuration of tlie Casimir energy measurement. Tlie radiation cavity 
bounded by two parallel perfectly-conducting plates separated by 21. The plate size is 
2L X 2L. 
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1 Introduction 

In the dawn of the quantum theory, the divergence problem of the specific heat of the radi- 
ation cavity was the biggest one (the problem of the blackbody radiation). It is historically 
so famous that the difficulty was solved by Planck's idea that the energy is quantized. In 
other words, the phase space of the photon field dynamics is not continuous but has the 
"cell" or "lattice" structure with the unit area {Ax- Ap) of the size 2TTh (Planck constant). 
The radiation energy is composed of two parts, Ecas and Ef^: 



where the parameter (3 is the inverse temperature, / is the separation length between two 
perfectly-conducting plates, and L is the IR regularization parameter of the plate-size. 
See FigHJ The second part Ejj is, essentially, Planck's radiation formula. The first one 
Ecas is the vaccuum energy of the radiation field, that is, the Casimir energy. It is a very 
delicate quantity. The quantity is formally divergent, hence it must be defined with careful 
regularization. (See App.B). Ecas/ {2LY does depend only on the houdary parameter /. 



EidEM — Ecas + Ef^ 





(1) 
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Figure 2: Graph of Planck's radiation formula. V{f3,k) = ^ k'^ / {e^'' ~ 1) (1 ^ 
(3 < 2, 0.01 < k < 10), A; (photon energy) = hw = hu = hc^, (3 = l/keT where 
ch = 2000 eV A, /cB(Boltzmann's constant) = 8.6 x 10~^ eV/K. k={l,10) eV correspond 
to ^ = (2000,200)A, ^=(1,2) eV-^ correspond to T=(10V8.6, 10Vl7.2) K. P = 0.1 corre- 
sponds to VAk = (2000)-^ X 0.1 X Ak[eV/k% 
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The quantity is a quantum effect and , at the same time, depends on the global (macro) 
parameter /. 

'""'^ ^ i?4(the fourth Bernoulli number) = , (2) 



(2L)2 (2/)3 4! 720 (2/)3 ' '30 

In Figj2], Planck's radiation spectrum distribution is shown. Introducing the axis of 
the inverse temperature (/3), besides the photon energy or frequency (k), it is shown stere- 
ographically. Although we will examine the 5D version of the zero-point part (the Casimir 
energy), the calculated quantities in this paper are much more related to this Planck's 
formula. El We see, near the /3-axis, a sharply-rising surface, which is the Rayleigh- Jeans 
region (the energy density is proportional to the square of the photon frequency). The 
damping region in high k is the Wien's region, [f] The ridge (the line of peaks at each (3) 



^ See the recent reviews on Casimir energy: Ref.[T]. 

^ Planck's formula depends only on the temperature not on I. (Comparatively the Casimir energy 
part does not depend on (3, but on the separation I. ) It is known that /3 can be regarded as the periodicity 
for the axis of the inverse temperature (Euclidean time). The axis corresponds to the extra axis in the 
following text. 

^ We recall that the old problem of the divergent specific heat was solved by the Wien's formula. This 
fact strongly supports the present idea of introducing the weight function (see SecIHI . 
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forms the hyperbolic curve (Wien's displacement law). When we will, in this paper, deal 
with the energy distribution over the 4D momentum and the extra-coordinate, we will see 
the similar behavior (although top and bottom appear in the opposite way). In order to 
compare this 4D case with the 5D case of the present paper, we do some preparation for 
the Casimir energy in App.B. 

In the quest for the unified theory, the higher dimensional (HD) approach is a fascinating 
one from the geometrical point. Historically the initial successful one is the Kaluza-Klein 
model[21 E], which unifies the photon, graviton and dilaton from the 5D space-time ap- 
proach. The HD theories , however, generally have the serious defect as the quantum 
field theory(QFT) : un-renormalizability. |fl The HD quantum field theories, at present, 
are not defined within the QFT. One can take the standpoint that the more fundamental 
formulation, such as the string theory and D-brane theory, can solve the problem. In the 
present paper, we have the new standpoint that the HD theories should be defined by them- 
selves within the QFT. In order to escape the dimension requirement D=10 or 26 from the 
quantum consistency (anomaly cancellation) we treat the gravitational (metric) field 
only as the background one. This does not mean the space-time is not quantized. See 
later discussions (SecI?]). We present a way to define 5D quantum field theory through the 
analysis of the Casimir energy of 5D electromagnetism. 

In 1983, the Casimir energy in the Kaluza-Klein theory was calculated by Appelquist 
and ChodosjS]. They took the cut-off (A) regularization and found the quintic (A^) di- 
vergence and the finite term. The divergent term shows the unrenormalizahility of the 
5D theory, but the finite term looks meaning full and, in fact, is widely regarded as the 
right vacuum energy which shows contraction of the extra axis. In this decade, triggered 
by the development of the string and D-brane theories, new treatments or new ideas were 
introduced to calculate the vacuum energy or the effective potential in HD. (The motiva- 
tion is to settle the stability problem of the moduli parameters jTJ |8]. ) One is to regard 
the system as the bulk and boundary, and do renormalization in both partsj9], [10]. Vari- 
ous regularization methods were carefully re-examined for the bulk-boundary theory [TT]. 
They are applied to various theories including realistic mo dels [T^ [TH| [HI [TSl US]. From 
the regularization viewpoint, the zeta-function (or dimensional) regularization combined 
with some summation formula is most commonly taken. The renormalization procedure, 
however, does not seem satisfactory. They succeed in calculating the properly regularized 
quantity and in separating the divergent terms. They found the finite part, but its physical 
meaning is obscure because the treatment of the divergent part is not established. They 
simply say, based on the analogy to the case of the ordinary (4D) renormalizable thoeries. 

The HD theories we consider here, are the HD generahzation of the famihar renormahzable (in 4D) 
theories such as 5D free scalar theory, 5D QED and 5D Yang-MiUs theory. 

^ Note that the ordinary power counting criterion is about the divergence degree for the couphng 
expansion. In this sense, there are some renormahzable HD field theories such as 6D $'^-theory. In the 
present paper, however, we have focused on the coupling-independent part, the Casimir part. This part is 
generally divergent in the higher dimensions. 

^ The gauge independence was confirmed in Ref.|S]. 
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the local counterterms can cancel divergences. lI They try to absorb divergences by the 
renormalization of parameters such as the brane tension (cosmological constant) and the 
gravitational constant. But it is fair to say that the A^-divergence problem, posed by Ap- 
pelquist and Chodos, is not yet solved. All these come from the unsatisfactory situation 
of the quantum treatment of the brane dynamics and the HD quantum field theories. 

In the development of the string and D-brane theories, a new approach to the renor- 
malization group was found. It is called holographic renormalization |T7 | IT8 | [T9 | |20 | |21 ] |22] . 
We regard the renormalization flow as a curve in the bulk (HD space). The flow goes along 
the extra axis. The curve is derived as a dynamical equation such as Hamilton- Jacobi 
equation. It originated from the AdS/CFT correspondence |23l \2M [25] . Spiritually the 
present basic idea overlaps with this approach. The characteristic points of this paper are: 
a) We do not rely on the 5D supergravity; b) We do not quantize the gravitational (metric) 
field; c) The divergence problem is solved by reducing the degree of freedom of the system, 
where we require, not higher symmetries, but some restriction based on the minimal area 
principle; d) No local counterterms are necessary. 

In the previous paper [26], we investigated the 5D electromagnetism in the flat geometry. 
For the later use of comparison (with the present warped case), we list the main results 
here. |^ The extra space is periodic (periodicity 21) and Z2-parity is taken into account: 

ds'^ = Tj^ydx^dx" + dy^ , — oo < x^, y < oo , y ^ y + 21, y ^ —y , 
= diag(-l, 1, 1, 1) , (X*0 = (X^ = X' = y)^ (x, y) , 

M,N = 0, 1, 2, 3, 5; fi,u = 0, 1, 2, 3. (3) 

The IR-regularized geometry of this 5D fiat space (-time) is depicted in Fig|21 The Casimir 
energy is rigorously (all KK-modes are taken into account) expressed as 

27r^ 

Ecas{I^,l) = -— dp dyp^W{p,y)F{p,y) , 

[^'^r Ji/i Ji/A 

N ^ /- X ,^4_/ N ~ — 3 cosh H2v — /) — 5 cosh A;/ 

F(p, y) = F"(p, y) + 4F+(p, y) = / dk \ / . (4) 

jp Ismhykl) 

where A is the 4D-momentum cutoff, and W{p, y) is the weight function to suppress the IR 
and UV divergences. |^ We obtained the following A and / dependence by the numerical 
analysis. 

^ See the first footnote of Sec.8 about the present treatment of the local counterterms. 

^ We do not require the reader to read ref.[26|. The necessary key procedures are explained in the text. 

^ Z2-odd part F~ comes from the quantum fluctuation of the extra component Ac^{X), while Z2-even 
part from the 4D components A^{X). 

The expression Q of Ecas is negative definite. The same thing can be said about the warped case 
(23) in the later description. 
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Figure 3: IR-regularized geometry of 5D flat space The 4D world (3-brane) is Eu- 
clideanized and is shown as 4D ball (shaded disk region) surrounded by sphere with 
radius fi~^. fi is the 4D IR regularization parameter, and is taken to be /i = l/l. UV- 
regularization is introduced, in SecJH by replacing the 4D ball with the "sphere lattice" 
composed of many small (size: 1/A) 4D balls. See Sec|5]for detail. 

t 5^" 
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1) Un-weighted case: W = 1 



Un-restricted integral region : 
EcasiA, = ^ r-0.1249/A^ - (1.41, 0.706, 0.353) x 10"^ lA^ ln(/A)l , 

Randall-Schwartz integral region : E/Jf^ = -^[-0.0894 A^] . (5) 

The quintic divergence of the upper one of shows the unrenormalizability of the 5D 
theory in the ordinary treatment. The triplet data show the unstable situation of numerical 
results. Ell As for the lower case, the (p, |/)-integral region is restricted to below the 
hyperbolic curve py = 
2) Weighted case 



-2.50^ + (-0.142, 1.09, 1.13) ■ 10"^^ for W = J-e-(i/2)«V-{i/2bV;^ = 
-6.04 10-2^ - (24.7, 2.79, 1.60) ■ 10"^^ for W = ^e"?^^ = W2 (6) 

-2.51^ + (19.5, 11.6, 6.68) ■ 10"^^ for W = ^e-(''/2)(P^+i/s/^) = 

where some representative cases {Wi: elliptic, iVi = 1.557/87r2 ; W2: hyperbohc, N2 = 
2(/A)^/87r^; Ws'. reciprocal, Ng = 0.3800/87r^) are shown. (See ref.|26] for other cases. ) 
The quantity Al is the normalization factor in the numerical analysis. The Casimir energy 
behavior of the case W2 is consistent with the Randall-Schwartz's one of 1). These results 
imply the renormalization of the compactification size I. 

<a./A/ = -^(l-4cln(/A)) = -^ ' ^=9(1^^4 = " ' 

where a and c should be uniquely fixed by clarifying the meaning of the weight function 
W and the unstable situation of the triplet data. 

The aim of this paper is to examine how the above results change for the 5D warped 
geometry case. The IR-regularized geometry of the 5D warped space(-time) is depicted in 
FigJH One additional massive parameter, that is, the warp (bulk curvature) parameter a; 
appears. The limit a; — )■ leads to the flat case. This introduction of the "thickness" l/u 
comes from the expectation that it softens the UV-singularity, which is the same situation 
as in the string theory. See Ref.|27j and [28] besides this work. 

The content is organized as follows. We start with the familiar approach to the 5D 
warped system: the Kaluza-Klein expansion, in Secj2l In Sec|3l the same content of Sec|2] 

" The results of (jH are based on the numerical integral of (g]) for ^ = (10, 20,40), A = 10 10^. The 
triplet coefficients correspond to the three values of I. This unstable situation does not appear in the 
present case of warped geometry. See ([5^ and (|116p . The same thing can be said about the weighted case 
2) in the following. 

This restriction was taken in Ref.[53 to suppress the UV-divergence. 
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Figure 4: IR-regularized geometry of 5D warped space ( ITOl) . The 4D world (3-brane) is 
Euclideanized and is shown as a 4D ball (shaded disk region) surrounded by sphere 
with radius fi~^. fi is the 4D IR regularization parameter. UV-regularization is introduced 
by replacing the 4D ball with the "sphere lattice" composed of many small (size: 1/A) 4D 
balls. See SecjSlfor detail. 
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is dealt in the heat-kernel method and the Casimir energy is expressed in a closed form 
in terms of the P/M propagator. The closed expression of Casimir energy enables us 
numerically evaluate the quantity in SecJH Here we introduce UV and IR regularization 
parameters in (4D momentum, extra coordinate)-space. A new idea about the UV and IR 
regularization is presented in SecJSl The minimal area principle is introduced. The sphere 
lattice and the renormalization flow are explained. In Secj6] an improved regularization 
procedure is presented where a weight function is introduced. Here again the minimal 
surface principle is taken. The meaning of the weight function is given in Sec|71 In Sec|H] 
we make the concluding remarks. Renormalization of the warp parameter u is explicitly 
shown. We argue the (4D) coordinates or momenta look quantized in the present treatment. 
The cosmological constant is addressed. We prepare five appendices to supplement the 
text. App.A deals with the classification of all minimal surface curves in the 5D warped 
space. App.B reviews 4D Casimir energy (the ordinary radiation cavity problem) where 
the features of the cut-off and zeta-function regularizations are examined. App.C explains 
the numerical confirmation of the (approximate) equality of the minimal surface curve and 
the dominant path in the Casimir energy calculation. The results, appearing in this paper, 
heavily relies on some numerical calculations. We explain them in App.D. Normalization 
constants of various weight functions are explained in App.E. 

2 Kaluza-Klein expansion approach 

In order to analyze the 5D EM-theory, we start with 5D massive vector theory. 

S5dv = j d'xdzy^i-^FMNF'''' - ^m'A^'AM) , Fmn = ^M^iv - O^Am , 

ds^ = -^h.^dx^'dx" + dz^) = GMNdX^^dX^ , G = det Gab , 
u^z^ 

(X^^) = (x^,^) , M,A^ = 0,1, 2,3, 5(or ^); = 0,1,2,3 . (8) 

The 5D vector mass, m, is regarded as a IR-regularization parameter. In the limit, m = 0, 
the above one has the 5D local-gauge symmetry. Casimir energy is given by some integral 

where the (modified) Bessel functions, with the index u = + appear. (See, for 
example, ref.|29j.) Hence the 5D EM limit is given by = 1 (m = 0). We consider, 
however, the imaginary mass case m = iu (m^ = — w^, z/ = 0) for the following reasons: 1) 
the UV-behavior does not depend on the bulk mass parameter m which is regarded as a IR 
regularization one; 2) we can compare the result with the 5D fiat case where the 5D scalars 
(4 even-parity modes + 1 odd-parity mode) are considered [26|, 3) u = Bessel functions 
are meaningfully simple in the analysis. Efl We can simplify the model furthermore. Instead 
of analyzing the = — of the massive vector ([8]), we take the 5D massive sca/ar theory 

^ 1/2 {iv? = — 3a;^/4) is another simple case where Bessel functions reduce to trigonometric 
functions. 
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on AdSs with = —Aoo^, v = + ni^/u"^ = 0. 

£ = v^(-Jv^$Va/$- Jm2$2) , G = detG 



MN 



V*^Va/$ - + J = 



(9) 



where ^{X) = ^{x'^,z) is the 5D scalar field. The background geometry is AdSs which 
takes the following form, in terms of z, 



MN 







1 



1 111 

— < z < or — < z < — 

T oj uj T 



{u;\z\y 



CO 

T- u 



(10) 



where we take into account Z2 symmetry: z ^ —z. uj is the bulk curvature (AdSs param- 
eter) and is the size of the extra space (Infrared parameter). 

In this section, we do the standard analysis of the warped system, that is, the Kaluza- 
Klein expansion approach. 

The Casimir energy Ecas is given by 



exp{i / d^XC} 



P<l>(X)exp 



I I d^xdz- ^ ^^r,,2.2;, ^a^ , ,un5 ^ 1 



{u\z\)^2 



^{u^z^dad^^ + iulzlY 



dz (uzY 



9^$ - m^^} 



.(11) 



Here we introduce, instead of $(X), the partially (4D world only) Fourier transformed 



field ^p{z). 



$(X) 



(2vr)^ 



(12) 



Eq. (fTT]) can be rewritten as 



V%{z) X 



exp 



d'^p 



2 dz^-^,iz){- 



1 2 d 1 d 



[UJZ 



IP + 



m 



dz {ojzY dz {(jjz) 



(13) 



where we have used the Z2-property, defined in ()15p . of ^p{z). From the above expression, 
we can read the measure function s{z) and the extra-space kinetic operator L^. 



d 



d 



m 



s{z 



[UZ 



dz {uzy dz (uz)^ ' 



(14) 
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and consider the Bessel eigen-value problem. 



Ipniz) = --Ipni-Z) for P 



{S{z)-'L, + Mn^}ijn{z)=0 
tpn{z) = tpn{~z) for P = + 



(15) 



with the appropriate b.c. at fixed points. Because the set {ipniz)} constitute the orthonor- 
mal and complete system, we can express ^p{z) as 



*p(^) = ^Cn{p)tpn{z 

n 

Eq. f|T3l) can be further rewritten as 
Vf^pi^z) exp 



(16) 



^ I (0i2 dz <j ^%{z)s{z){s{zr'L, - p')%{z] 



exp 



^ ' n 

exp$^{-^ln(p| + M„^)} ,(17) 



n,p 



where the orthonormal relation 



tpn{z)s{z)lljmiz)dz = S„ 



is used. This shows that s{z), defined in (1141) . plays the role of "inner product measure" 
in the function space (2;), l/w < z < l/T}. In (ITTl) . Wick's rotation is done for the 
time-component of {p'^}. 

• n 4 

tp ^ p , 



idp dp dp dp — )■ dp dp dp dp = d pe 



P 



The expression f[T7j) is the familiar one of the Casimir energy. 
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3 Heat-Kernel Approach and Position/Momentum Prop- 
agator 

Eq. (fT71) is the expression of Ecas by the KK-expansion. In this section, the same quantity 
is re-expressed in a closed form using the heat-kernel method and the P/M propagator. 



T)^p{z) exp 



l/T 
1/ui 



exp 



exp 



T- 



rp — 3 

d^PE 
(2vr)4 



dz{-<!>,{z)s{z){s{zr'L,-p')<^,{z] 



(2vr)4- 

l/T 



dzs{z) 



1/lu 



:gt(s(^)-iL,-p|)^^ + const 



(20) 



where we have used a formula 



dt = In M , det M > , M : a matrix 



(21) 



(The factors T"^ and T'^ in ([20]) come from the dimensional analysis. T has the meaning 
of the renormalization point. ) The above formal result can be precisely defined using the 
heat equation. 



e ^ ^^"^ = (const) X exp 



T' 



d^PE_ 
(27r)4 



~ Idt 
27 



Tr Hpj^{z,z';t) 



Ti Hpiz,z';t)= siz)Hp{z,z;t)dz , {-- (s^'L, - p^)}Hp{z, z';t) = . (22) 

The heat kernel Hp[z, z'\ t) is formally solved, using the Dirac's bra and ket vectors {z\, \z), 
as 



Hp{z, z'] t) = {z\e 



-(-s-iL^+p2)t| / 



(23) 
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(The bra and ket vectors {z\, \ z) are precisely defined by the orthonormal and complete set 
of Lz'. {ipniz)}- ) Using the set {'ipn(z)} defined in (fT5|) . the explicit solution of (122!) is 
given by 

neZ 

E,{z,z';t) = Y,^-^'''-^^'^'l{M^)M^')+^^^^^^^ , P = + , (25) 

nGZ 

where we have used the dimensionality of Hp and Ep read from fl22l) . {[Ep\ = [Hp\=L~^) . 
The above heat-kernels satisfy the following b.c. 



nGZ 

\im^Ep{z,z';t) = J2l{Mz)Mz') + Mz)M-z')} = i^\z\?S{\z\ - \z'\) , P = + 



neZ 



where e{z) is the sign function. The above equation defines (5(|2| — We here introduce 
the position/momentum propagators as follows. 

G;{z,z')= / dt Hp{z,z';t) = ... -{ij,,{z)Mz')-Mz)M-^')} , 

poo 

G;{z,z')= dt Ep{z, z'; = ^ , , -{Mz)Mz') + Mz)M-^')} ■ (27) 
^0 nez^*'^"+^ ^ 



1 /I 



(1) Definition il\i{z) — {n\z) ~ (-^1^^) , 
(2) Z2 - property | - z) P\z) , (-z| = P(z| , P - t1 

(3) Orthogonality 

dz dz 

{n\z){z\k) ^2 —{n\z){z\k) = 5n,k, {n\k) = 4,fc , 



(^1^1)3^ ' ' ' 7i (loz) 
{z\z') 



{uj\z\)h{z)e{z')6{\z\-\z'\) for P= - 1 
{u;\z\f5i\z\-\z'\) for P=l 

(4) Completeness 



1 /I 



;l-)(-l = 2 /"-^|z)(z| = l 



(tj|z|)3' ' Jj_ iujz)^' 

Y,\n)in\ = l , (24) 
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They satisfy the following differential equations of propagators. 

(4 - p's{z))G;iz, z') = J2 T 

^ r eiz)e{z')6i\z\ - \z'\) for P= - 1 

\ S{\z\ - \z'\) for P=l ^ ' 

Therefore the Casimir energy Ecas is, from (122|) and f l25|) . given by 

where s(2;) = 1/{ujzY (IT^ . This expression leads to the same treatment as the previous 
section. Note that the above expression shows the negative definiteness of E^^^. 
Here we introduce the generalized P/M propagators, Ia{P=—) and Jq(P=+) as 

^a{p'^;z,z') = —Hp{z,z';t) 

Ja{p'^;z,z') = / —Ep{z,z';t) 
Jo 

r°° r/f 1 

= / t;;E^"^*'""'''^*9{^-W^"(^')+^"W^"(-^')> ' P=+ ' (30) 

where a is the arbitrary real number. Then we have the following relations. 

Io{p'^;z,z') = G~{z,z') , Jo{p'^;z,z') = G'^{z,z') 

^^"'^^ = -Ia-l{p'^;Z,z') , [ dk'^Ia{k'^]Z,z') = Ia+l{p'^]Z,z') 

dJa{p ,Z,Z)_ ^ , [ dk'^Ja{k'^]Z,z') = Ja+l{p'^]Z,z') 



dp 

I 

1 ' 

dp"^ Jp^ 
[p^ - r{z)-^L^)Ip{p'^; z, z) = - (3 1 p+iip"^ ; z, z) , (p^ - r{z)~^L^)Jfi{p^; z, z') = -/3J/3+i(p^; z, z) 

dUp-'-z.z') r^^t 

— ' -dtHplz, z ; t) 



da L t 



We notice the subtraction of positive infinity (M-independent term) in the formula (PT|) is essential 
for this negative definiteness. This should be compared with the expansion-expression Ecas of (|17p . 
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Finally we obtain the following useful expression of the Casimir energy for P = =]=. 



d^PE 
(27r)4 



d^PE 



{TTGt{z,z')}dk^ 



d^PE 



i/r ^oo 
dz s(z 



vl 



(2vr)4 7,1 

Here we list the dimensions of various quantities appeared above 



{Yi.J^{k^;z,z)}dk^ 
{Gl{z,z)\dk^ . 



(32) 















L 








-l-ai J a 


Ecas 


$ 


A,p,uj,T 
H E 

6{z - z') 




s{z) 
\n), {n\ 


z, I 

y-j p 


t 




Cn{p) 



(A is a regularization parameter defined below.) 

The P/M propagators GJ in (IHT]) and f l32p can be expressed in a closed form. 

(See, for example, jS].) Taking the Dirichlet condition at all fixed points, the expression 
for the fundamental region {l/u < z < z' < 1/T) is given by 

u;3^2.,2{Io(!)Ko(pz) tKo(5)Io(;5^)}{Io(|)Ko(p^') TKo(|)Io(;^^')} 



G^{z,z') = T^z-z 



Io(|)Ko(5)-Ko(|)Io(f; 



P 



p^ > (space-like) .(33) 



We can express the A-regularized Casimir energy in terms of the following functions 
F^ip,z). 

-l/T 



EiZ{uj,T)= I I dzF^{p,z 



[UZ 



(2^)1p<A^l/a 

F^{p,z)^s{z) / {GUz,z)}de 
kG]:{z,z)dk= / J^^{k,z)dk , 



(34) 



where T^{k,z) are the integrands of F^{p,z) and p = \[p^- Here we introduce the UV 
cut-off parameter A for the 4D momentum space. In Figj5]and FigEl we show the behavior 
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Figure 5: Behavior of In 2)] = ln|^ G]^{z,z)/ {uzf\. u = 10^T = 1,A = 2 x 10^ 

1.0001/a; < z < 0.9999/T. AT/co < k < A. Note In |(l/2) x (1/2)| ^ -1.39. 




of J-'^(k,z). The table-shape graphs say the "Ray ley- Jeans" dominance. That is, for 
the wide-range region (p, z) satisfying both p{z — ^) ^ 1 and — 2;) ^ 1, 

{p,z) e{{p,z)\p{z--):^i andpi^- z):^l} . (35) 



4 UV and IR Regularization Parameters and Evalua- 
tion of Casimir Energy 

The integral region of the above equation (IMl) is displayed in FigJTl In the figure, we 
introduce the regularization cut-offs for the 4D-momentum integral, /i < p < A. As for the 
extra-coordinate integral, it is the finite interval, l/u < z <1/T = e'^''/u, hence we need 
not introduce further regularization parameters. For simplicity, we take the following IR 

The energy density in (p", z)-space is approximately given by, using F^{p,z) w —^p+ ^A. For 
small p, F^{p,z) sa ^ A (const.). This should be compared with Ep of (fT07| : ^/(c'^" - 1) - 1//3 (const) 
for small cj. 



16 



Figure 6: Behavior of In 2)] = ln|^ Gl{z,z)/ {ujzf\. u = 10^T = 1,A = 2 x 10^ 

1.0001/a; < z < 0.9999/T. AT/co <k<A 




cutoff of 4D momentum. : 

T 



^ = A . _ = Ae— " . (37) 



Hence the new regularization parameter is A only. 
Let us evaluate the (A, T)-regularized value of (134|) . 



Ec'Zi^,T) = ^^ j\pj'\zp'F^{p,z) 



F^iP,z) = j^^l_^~kG^{z,z)dk . (38) 

The integral region of {p, z) is the rectangle shown in FigJT] . 

Note that eg. (138!) is the rigorous expression of the (A, T)-regularized Casimir energy. We 
show the behavior of {—l/2)p^F~{p, z) taking the values u = 10^, T = 1 in FigIH](A = 10^), 

If we take the following relation furthermore 

A = cj , (36) 

then fi ~ Tand we need not any additional regularization parameters. We do not take this relation. The 
choice of the regularization parameters affects the counting of the divergence degree. See later discussion 
of ea.(|il|) 
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Figure 8: Behaviour of {-l/2)p^F-{p,z) 038]). T = l,w = 10^A = lOl l.OOOl/w <z< 
0.9999/T, AT/cj < p < A . 
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Figure 9: Behavior of {-l/2)p^F"{p, z) T = 1,lo = 10^A = 2 ■ 10^ 1.0001/u; < 

z < 0.9999/T, AT/u <p<A. 




Figure 10: Behavior of {-l/2)p^F-{p, z) T = l,u = 10^A = 4 ■ lOl 1.0001/w < 

z < 0.9999/T, AT/u <p<A. 
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Figini^A = 2 ■ 10'') and FigJTOlfA = 4 ■ 10^). All three graphs have a common shape. 
(We confirm the graphs do not depend on the choice of u and T very much.) Behavior 
along p-axis does not so much depend on z. A valley runs parallel to the 2;- axis with the 
bottom line at the fixed ratio of p/A ~ 0.75. The depth of the valley is proportional 
to A^. Because Ecas is the (p, z) 'fiat-plane' integral of p^F{p, z) , the volume inside the 
valley is the quantity Ecas ■ Hence it is easy to see Ecas is proportional to A^. This is the 
same situation as the fiat case (the upper eq. of ([5])). Importantly, (138|) shows the scaling 
behavior for large values of A and 1/T. From a close numerical analysis of (p, 2;)-integral 
we have confirmed 



2t^ 
W) 



A' 

-0.0250 — 
T 



(39) 



which does not depend on u and has no In —term. (Note: 0.025 = 1/40. See App.D for the 
numerical derivation.) Compared with the flat case (the upper eq. of ([5])), we see the 
factor plays the role of IR parameter of the extra space. We note that the behavior of 
FigJSllTUjis similar to the Rayleigh- Jeans 's region (small momentum region) of the Planck's 
radiation formula (Figj2]) in the sense that p^F{p^ z) oc p^ for p <^ A. 

Finally we notice, from the FiglHlfTOl the approximate form of F{p^ z) for the large A 
and 1/T is given by 

F^(p,2)^^A(l-|) , / = 1 , (40) 

which does not depend on 2, uj and T. / is the degree of freedom. The above result is 
consistent with fl35l). 



5 UV and IR Regularization Surfaces, Principle of Min- 
imal Area and Renormalization Flow 

The advantage of the new approach is that the KK-expansion is replaced by the integral 
of the extra dimensional coordinate z and all expressions are written in the closed ( not 

The requirement for the three parameters oj, T, A is A a; T. See ref.[32] for the discussion about 
the hierarchy A,Ld,T. In the apphcation to the real world, the most interesting choice is T ITeV = 
103GeV(TeV physics), w - 10i^GeV(GUT scale), and A - 10i^GeV(Planck mass,Mp,). In the numerical 
calculation, however, we must be content with the appropriate numbers, shown in the text, due to the 
purely technical reason. Another interesting choice is w ~ 10^'^eV(neutrino mass, rrii, ^ y/ Mpi / Rcos ,Rcos'- 
cosmological size), A - lO^^GeV = lO^^eV (Planck mass) and T - lQ-'^^eY{r^ R^-}{MpiRcosY/^), See the 
discussion about the cosmological term in the concluding section. 

The Valley-bottom line 'path' p = p{y) k, 0.75A corresponds to the solution of the minimal principle: 
5Si = 0, Si[p{z),z] = (l/87r2) / dzp{zf F {p{z) , z) , F ^ -{A-p)/2 (gO]). This will be referred in SecEI. 

■^^ The result ([5^ is based on the numerical calculation for the following cases: 1) T = l,w = 10"', A = 
104x(l,2,4,8, 16); 2)T = l,a; = 10^ x (1, 2, 4, 8, 16), A = 2x10^3) T = (1, 2, 4, 8, 16), a; = 10^, A = 2x10''. 

This numerical result can be checked using the approximate form (PPI) . Ii/^ dz p^-{—}^){A—p) = 

-0.025^(1 + 0(r/a;)). 
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expanded ) form. The A^-divergence, (139|) . shows the notorious problem of the higher 
dimensional theories, as in the flat case (the upper eq. of ([5])). In spite of all efforts of the 
past literature, we have not succeeded in defining the higher-dimensional theories. (The 
divergence causes problems. The famous example is the divergent cosmological constant 
in the gravity-involving theories. j5] ) Here we notice that the divergence problem can be 
solved if we find a way to legitimately restrict the integral region in (p, z )-space. 

One proposal of this was presented by Randall and Schwartz |29]. They introduced 
the position- dependent cut-off, fi < p < A/uu , u G [l/o;, 1/T] , for the 4D-momentum 
integral in the "brane" located at z = u. See FigJTl The total integral region is the lower 
part of the hyperbolic curve p = A/ujz. They succeeded in obtaining the finite /3-function 
of the 5D warped vector model. We have confirmed that the value Ecas of (138|) . when 
the Randall-Schwartz integral region (Fig|7]) is taken, is proportional to A^. The close 
numerical analysis says 

E^^fiu,T) = — — dq dz q'F~iq,z) = -— / du / dpp^F^iP^u) 

-1.58 X 10"^ - 1.69 X 10"Mn - ^ ,(41) 



(27r)4 UJ { ' ' UJ 

which is independent of T . E§ Efl This shows the divergence situation does not improve 
compared with the non-restricted case of fl5^ . T of is replaced by the warp parameter 
UJ. This is contrasting with the flat case where E^^^ cx — A^. (the lower eq. of ([5])) The 
UV-behavior, however, does improve if we can choose the parameter A in the way: A oc w. 
This fact shows the parameter uj "smoothes" the UV-singularity to some extent. 

Although they claim the holography is behind the procedure, the legitimateness of 
the restriction looks less obvious. We have proposed an alternate approach and given a 
legitimate explanation within the 5D QFT|31l l33l l26l IM] . Here we closely examine the 
new regularization. 

On the "3-brane" at 2; = l/cj, we introduce the IR-cutoff fx = A - - and the UV-cutoff 
A (^ < A). See FigHH 

/i < A (T < w) . (42) 

This is legitimate in the sense that we generally do this procedure in the 4D renormalizable 
thoeries. (Here we are considering those 5D theories that are renormalizable in "3-branes". 



^•^ The approximate form (^0)) predicts the similar result. dq /j^^"' dz ■ (— i)(A — (7) = ~fj^ ^(1 + 
0((T/w)3)). 0.01666 • • • = 1/60. 

23 The result (gl]) is based on the numerical-integral data for T = (1, 2, 4, 8, 16), a; 10^, A = 2 x lO''; 
A = 10^ X (l,2,4,8,16),r = 1,0; = 10^; cj = 10^ x (1, 2, 4, 8, 16), T = 1,A = 2 x lO'^. See App.D for the 
numerical derivation. 

■2'* In ref.[29], they take A = O.Sw, w, 2aj, • • • and evaluate /3-function (of the gauge coupling constant) for 
the different cases. They regard the parameter ui as the physical cutoff. This choice, however, is not allowed 
in the present standpoint A 3> w ^ T. The fact that u appears as (j4ip imply the warp parameter can 
control the UV-behavior to some extent. It matches the belief that the theoretical parameter oj physically 
means the extendedness of the system configuration and smoothes the UV-singularity. 
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Figure 11: Space of (p,z) for the integration (present proposal). 




Examples are 5D free theories (present model), 5D electromagnetism|26| . 5D $^-theory, 
5D Yang-Mills theory, e.t.c.) In the same reason, on the "3-brane" at z = 1/T, we may 
have another set of IR and UV-cutoffs, /i' and A'. We consider the casj^: 

/i' < A', A' < A, /i ~ /i' . (43) 

This case will lead us to introduce the renormalization flow. (See the later discussion.) 
We claim here, as for the regularization treatment of the "3-brane" located at other points 
z {l/bJ < z < 1/T), the regularization parameters are determined by the minimal area 
principle. @ To explain it, we move to the 5D coordinate space {x^^z). See FigJT2l The 
p-expression can be replaced by ^x^^-expression by the reciprocal relation. 

^Jx^{z)x^'{z) = r{z) ^ ^ . (44) 

The UV and IR cutoffs change their values along ^-axis and their trajectories make surfaces 
in the 5D bulk space (x'^jz). We require the two surfaces do not cross for the purpose 
of the renormalization group interpretation (discussed later). We call them UV and IR 

Another interesting case is /i < A, A <C A', /i ^ fi' . This case gives us the opposite direction flow. 

We do not quantize the (bulk) geometry, but treat it as the background. The (bulk) geometry 
fixes the behavior of the regularization parameters in the field quantization. The geometry infiuences the 
"boundary" of the field-quantization procedure. 
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Figure 12: Regularization Surface Bm and Buv in the 5D coordinate space The 
three graphs at the bottom show the flow of coarse graining (renormahzation) and the 
sphere lattice regularization which will be explained after some paragraphs. 



X 




r„„ (z) 



regularization (or boundary) surfaces(5(/y, 5/^). 

Buv ■■ V(xi)2 + (x2)2 + (x3)2 + (a;4)2 = ruv{z) , -<z<l/T , 

UJ 

BiR : + (x2)2 + (x3)2 + (x4)2 = riR{z) , -<;s<l/T , (45) 

UJ 

where ruy{z) and rjR{z) are some functions of 2 which are fixed by the minimal area 
principle. The cross sections of the regularization surfaces at z are the spheres with the 
radii ruy[z) and rjR{z). Here we consider the Euclidean space for simplicity. The UV- 
surface is stereographically shown in FigHSjand reminds us of the closed string propagation. 
Note that the boundary surface Buv (and Bir) is the 4 dimensional manifold. 

The 5D volume region bounded by Buv and Bjr is the integral region of the Casimir 
energy Ecas- The forms of ruv{z) and rjR{z) can be determined by the minimal area 
principle. 

^(Surface Area) = , 

3+-rV-^— = , r' = ^ , r" = ^ , l/uj<z<l/T . (46) 
z r'2 + 1 dz ' dz^ ' - - ' ^ ' 

In App.A, we present the classification of all solutions (paths). It helps to find appropriate 
minimal surface curves for the renormahzation flow. 
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Figure 13: UV reguiarization surface (Buy) in 5D coordinate space. 




£=1/(jO 



radiu s 1 /A' 



UV,1/T 



In FigJH] we show two result curves of 
(r' = dr/dz) : 



taking the following boundary conditions 



Fig.ll: Coarse Conf. goes to Fine Conf. as z increases 
IR-curve (upper): r[l] = 0.8,r'[l] = 1.0 type (ia) 
UV-curve (lower): r[l] = 10~^,r'[l] = —1.0 type (ia) 



(47) 



They show the flow of renormalization really occurs by the minimal area principle. 
(See the next paragraph for the renormalization flow interpretation.) These results imply 
the boundary conditions determine the property of the renormalization flow. 

The present reguiarization scheme gives the renormalization group interpretation to 
the change of physical quantities along the extra axis. See FigJT^ PI In the "3-brane" 
located at z, the UV-cutoff is ruv{z) and the reguiarization surface is the sphere 5*^ with 
the radius ruv{z). The IR-cutoff is rm^z) and the reguiarization surface is the another 
sphere 5''^ with the radius rm^z). We can regard the reguiarization integral region as the 



The flow direction is opposite to the one shown in FigfTTl 

The minimal area equation (j46p is the 2nd derivative differential equation. Hence, for given two 
initial conditions (;for example, r{z = 1/w) and dr / dz\z=i/S}^ there exists a unique solution (path). The 
presented graphs are those with these initial conditions. Another way of choosing the initial conditions, 
r{z ~ 1/w) and r{z = 1/T), is possible. Generally the solution of the second derivative differential equation 
is fixed by two initial conditions. 

This part is contrasting with AdS/CFT approach where the renormalization flow comes from the 
Einstein equation of 5D supergravity. 
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Figure 14: Numerical Solution by Runge-Kutta. (H6|l . Vertical axis: r; Horizontal axis: 
z. T = l,u = 10^10-^ < z < 1.0. Upper (Bm): r(l) = 0.8,r'(l) = 1.0,; Lower (Buv)- 
r(l) = 10~^,r'(l) = —1.0 . Both curves are Graph Type (ia). 




0.2 0.4 0.6 0.8 1 



sphere lattice of the following properties: 

A unit lattice (cell) : the sphere with radius ruv{z)8ind its inside , 
Total lattice : the sphere S'^with radius rm^z) and its inside . 

It is made of many cells above , 

TT I • (48) 

ruv{z)J 

The total number of cells changes from (^)^ at z = 1/uj to (^)^ at z = 1/T. Along the 
z-axis, the number increases or decreases as 

''^^^'^^'^Niz) . (49) 



For the "scale" change z ^ z + Az, N changes as 

A(lnAr)=4|-{ln(^i44)}-^^ • (50) 
dz ruv{z) 

When the system has some coupling g{z), the renormalization group /3((7)-function (along 
the extra axis) is expressed as 

A(ln^) _ 1 _ 1 1 Idg 

^ A(\nN) A(lniV) g 4 IlMi^) g dz ' ^ ^ 
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where g{z) is a renormaUzed coupUng at z. r°l 

We have explained, in this section, that the minimal area principle determines the flow 
of the regular izat ion surfaces. 



6 Weight Function and Casimir Energy Evaluation 

In the expression (132|) . the Casimir energy is written by the integral in the {p, z) -space 
over the range: l/oo < z < 1/T, < p < oo. In SecJS], we have seen the integral region 
should be properly restricted because the cut-off region in the 4D world changes along the 
extra-axis obeying the bulk (warped) geometry {minimal area principle). We can expect 
the singular behavior (UV divergences) reduces by the integral-region restriction, but the 
concrete evaluation along the proposed prescription is practically not easy. In this section, 
we consider an alternate approach which respects the minimal area principle and evaluate 
the Casimir energy. 

We introduce, instead of restricting the integral region, a weight function W{p, z) in the 
{p, 2;)-space for the purpose of suppressing UV and IR divergences of the Casimir Energy. 



u;,T) 



POO 

F^{P,z) = s{z) / {Gl{z,z)}dk 



Examples of W{p^ z) : W{p, z 

-l„-(l/2)p2/'^^-(l/2)^^T2 



(iViJ 

(Nit) 



-lp-(l/2)pV^' 



= Wi{p,z), Ni = 1.711/87r2 
WibiP,z), iVife = 2/87r2 



(iV4 



{N2)-'e-P'''/^ = W2{p,z), N2 = 2^/871^ 
^ ie-(V2)?^^^^TW^iy3(p,^), iV3 = |^/8: 



{N5 



-1, 



r'e-P/^'-''' ^W,{p,z), iV5 = i/87r2 
We{P,z), iV6 = i/8vr2 
WjiP,z), Nj = 1/871' 
E Wsip^z), Ns = 0.4177/87r2 



-lg-p2/22a;2T 
-lp-(l/2)pV^^ 



(iV47)- 
(iV56 



'lp-l/2{p'/i^'){r/i^'+l/z'T') 



W,-r{p,z), K 
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g-i/2{p/.c.r)(p/c.+i/.T) ^ w^eiP^z), N^e 

L„-l/2{p2/'^'+l/^'T2)2 



0.1028/87r2 
= 0.1779/87r2 
0.01567/87r2 



P 



(iV9)-ie-i/2(p/<^+i/.T)2 ^ Wg{p,z), Ng = 0.05320/87r2 



iuz 



pI + pI + pI+ pI 



k G],{z, z)dk 



elliptic suppr. 
kinetic-energy suppr. 
hyperbolic suppr. 1 
hyperbolic suppr. 2 

linear suppr. 
parabolic suppr. 1 
parabolic suppr. 2 
higher-der. suppr. 1 
reciprocal suppr. 1 
higher-der. suppr. 2 
reciprocal suppr. 2 
higher-der. reciprocal suppr. 
reciprocal suppr. 3 



where ^^(z, z) are deflned in ( l33l) . The normalization constants Ni are explained in App.E. 



Here we consider an interacting theory, such as 5D Yang-Mills theory and 5D theory, where the 
coupling g{z) is the renormalized one in the '3-brane' at z. 
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Figure 15: Behavior of {—Ni/2)p^Wi{p, z)F (p, z)(eniptic suppression). A = 20000, u = 
5000, T = 1 . LOOOl/cj <z< 0.9999/T, = AT/u < p < A. 




Ell In the above, we hst some examples expected for the weight function W{p, z). W2 and 
W3 are regarded to correspond to the regularization taken by Randall-Schwartz. How 
to specify the form of W is the subject of the next section. We show the shape of the 
energy integrand {—l/2)p^W{p, z)F~{p, z) in Fig JTSlfTSl for various choices of W. We notice 
the valley-bottom line p 0.75A, which appeared in the un- weighted case (FiglSl fTUj) . 
is replaced by new lines: p"^ + z"^ x oj^T^ ~ const (Fig JT^ VTi ). pz const (FigJIBlVTci). 
p ^ const X z(FigJT7liy^). p ~ const x ^/z (FigJTSlWfi). They are all located away from 
the original A-effected line {p ~ 0.75 A), o 

We can check the divergence (scaling) behavior of E'^Jf by numerically evaluating the 

■^^ In the warped geometry we have, besides the cut-off parameter A, two massive parameter T and cu. 
We make aU exponents in ([5^ dimensionless by use of T for z, and uj for p. 

■^^ For the graphical view, we take rather large values of w's. If we take a more smaller value for uj, the 
position of a valley-bottom line deviates more from that of the un- weighted case {p ^ 0.75A). 
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Figure 16: Behavior of {—N3/2)p^W3{p, z)F (p, 2;) (hyperbohc suppression2). A = 
20000, cu = 5000, T = 1 . l.OOOl/o; <z< 0.9999/T, ^ = AT/u < p < A. 




{p, 2;)-integral ( 152|) for the rectangle region of FigJZl B 



w 

Cas 



T 
^ „ 

Ka' 

T. 



A 

rw.wu 111- 
-0.01 In^ 
0.08 In^ - 
0.07 In^ - 
0.07 In'^^ - 



^A X 1.2(1 + 0.11 In^ 
^A X 2.0(1 + 0.07 In^ 
X 0.062 (1 + 0.03 In 
X 0.14(1 
A X 1.5 (1 4 
^A X 1.5 (H 

^A X 0.86 (1 + 0.07 In"^ - 0.07 In'^ j for We (53) 
^^A X 1.3(1 + 0.06 In I 
I^A X 1.6(1 + 0.09 In I 
^A X 0.47(1 + 0.05 In^ 
^A X 0.93 (1 + 0.06 In^ 
^A X 1.1 (1 + 0.06 In^ 
^A X 0.91 (1 + 0.05 In'^^ 

The data fitting is based on the numerical integration for the different cases of (A, w, T). For example, 
the Wi formula is based on the numerical values of Ecas for T = 0.01, w = 10^, A = 10^ x (1, 2, 4, 8, 16); a; = 
10^A = 2 X 104,r = (1,1/2,1/4, 1/8, 1/16); T 1,A ^ 8 x 10^w 10^ x (8,16,32,64,128). See App.D 
for the numerical derivation. 
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Figure 17: Behavior of (—N4/2)p^W4{p,z)F (p, 2)(linear suppression). A = 100, u = 
10, T = 1 . 1.0001/w <z< 0.9999/T, = AT/u < p < 25. In order to demonstrate the 
valley-bottom line is similar to a minimal surface line (See App.C), we here take rather 
small values of A and u. The contour of this graph will be shown later. 
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Figure 18: Behavior of {—Nq/2)p^Wq{p,z)F (p, z)(parabohc suppression2). A = 
20000, u = 5000, T = 1 . l.OOOl/cj <z< 0.9999/T, /x = AT/cu < p < A. 




The suppression behaviors of W2 and improve, compared with f HT]) by Randall- 
Schwartz. The quintic divergence of ( 14T|) reduces to the quartic divergence in the present 
approach of W2 and W3. The hyperbolic suppressions, however, are still insufficient for the 
renormalizability. After dividing by the normalization factor, AT~^, the cubic divergence 
remains. The desired cases are others. The Casimir energy for each case consists of three 
terms. The first terms give finite values after dividing by the overall normalization factor 
AT^^. The last two terms are proportional to log A and show the anomalous scaling. Their 
contributions are order of 10~^ to the first leading terms. The second ones (In ^) contribute 
positively while the third ones (In ^) negatively. 

They give, after normalizing the factor A/T, only the log-divergence. 

E^^JAT"^ = -aoj^ (1 - 4cln(A/u;) - 4c'ln(A/r)) , (54) 

where a, c and c' can be read from (153|) depending on the choice of W . This means the 
5D Casimir energy is finitely obtained by the ordinary renormalization of the warp factor 
LO. (See the final section.) In the above result of the warped case, the IR parameter / in 
the fiat result ([7]) is replaced by the inverse of the warp factor u. 

At present, we cannot discriminate which weight is the right one. Here we list char- 
acteristic features (advantageous (Yes) or disadvantageous(No), independent (I) or depen- 
dent (D), singular (S) or regular (R)) for each weight from the following points. 

We should note that, for so many different forms of the suppression factor, the normalized E^^^ takes 
this log-divergence behavior ((M)) with similar coefRcients (a, c, c'). Exceptions are W2 and W3. 
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point 1 The behavior of W for the limits: T — oo (4D hmit) ; T — (5D hmit) ; 

— ^ (fiat limit). This property is related to the continuation to the ordinary 
field-quantization. 

point 2 flow the path (bottom line of the valley) depends on the scales T and u. 

point 3 Regular(R) or singular(S) at 2; = 0. This point is not important because the 
range of 2; is l/oo < z < 1/T or —l/T<z<~l/uj. 

point 4 Symmetric for p/u Tz. 

point 5 Symmetric for p/u ^ 1/Tz. (Reciprocal symmetry) 

point 6 The value of a. 

point 7 The values of (-4c,-4c'). 

point 8 Under the Z2-parity z O —z, W{p,z) is even (E), odd (O) or none (N). 
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W type 
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We notice Wi and Wg are specially important. 

So far as the legitimate reason of the introduction of W{p, y) is not clear, we should 
regard this procedure as a regularization to define the higher dimensional theories. We give 
a clear definition of W{p,y) and a legitimate explanation in the next section. It should be 
done, in principle, in a consistent way with the bulk geometry and the gauge principle. 



7 Meaning of Weight Function and Quantum Fluctua- 
tion of Coordinates and Momenta 

In the previous work|26|. we have presented the following idea to define the weight function 
W{p,z). In the evaluation (l52l) : 



- E'^^X^^T) = dz W{p, z)F^p, z) 

= J^jdpJ^^ dzfW{p,z)F^{P,z) , (55) 

the (p, 2;)-integral is over the rectangle region shown in FigUl] (with A — ?■ oo and /i — )■ 0). 
F^[p, z) is explicitly given in fjM|) . Following Feynman|35], we can replace the integral by 
the summation over all possible pathes p{z) as schematically shown in FigJT9l 
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Figure 19: The mesh of the dotted hnes shows the the ordinary integral of J dpdz. Two 
sohd hnes show two pathes pi{z) and p2{z)- The path-integral / Vp{z) is the integral over 
all possible pathes. 
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S\p{z), z] = —p{zfW{p{z), z)F^{p{z), z) . (56) 

Especially, in the figure, the mesh shows the independency of the integral-variables and 
z. Two pathes pi{z) and P2{z) are shown as two solid lines. There exists the dominant 
path pw{z) which is determined by the minimal principle E§ : 55 = 0. 



dp 



d\n{WF) 



Dominant Path pw{z) : ^ = Q^^y^ ■ (57) 

p ~^ dp 

Hence it is fixed by W{p^z). Examples are the valley-bottom lines in Fig JTMTSl On the 
other hand, there exists another independent path: the minimal surface curve rg{z). 

4 r"r 1 1 

Minimal Surface Curve rJz) : 3 H — r'r ^ = , — < z < — , (58) 

z r' + 1 bJ T 



The valley-bottom line of the graph S = S[p, z] can be obtained by two steps. First we take the 
expression ([55]) . — i?^^(w,T) = J dp J^j^ dz S[p,z], S[p,z] oc p^W{p, z)F{p, z). We do the variation, 
assuming the two independent coordinates p and z: z — >■ z + Az, p — )■ p+Ap, S — ^ S + AS, AS — ApdpS + 
AzdzS. Secondly we put the condition of path: p = p{z). Then Ap = Az, AS* = i^dpS + dzS)Az 
From the variation condition AS* = 0, we obtain ([5 
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which is obtained by the minimal area principle: 



ds^ = „ „ (Sab + 7 r7:)dx"'dx'' = QabMdx^'dx'' 

irr')'^ 



5 A = , A = \ Vdet(?a6 d^x= \ — 

J Ji/lo ^ ' 



1 r^dz . 



(59) 



See App.A for detail. Hence rg{z) is fixed by the induced geometry Qabi^)- Here we put 
the requirement^ 



(60) 



where pg = l/vg. This means the following things. We require the dominant path coincides 
with the minimal surface line Pg{z) = l/rg{z) which is defined independently of W{p,z). 
In other words, W{p,z) is defined here by the induced geometry Qabi^)- In this way, we 
can connect the integral-measure over the 5D-space with the (bulk) geometry. We have 
confirmed the (approximate) coincidence by the numerical method. (See App.C) 

In order to most naturally accomplish the above requirement, we can go to a new 
step. Namely, we propose to replace the 5D space integral with the weight W, (155|) . by 
the following path-integral. We newly define the Casimir energy in the higher-dimensional 
theory as follows. 

-£casiu;,T,A) = dp [ T\Vp''iz)Fip,z) 

Jl/A Jp{l/u)=p{l/T)=l/p 



X exp 

1/m 



1 f^'^ 1 1 



1 dz 



dp! T\Vx'^{z)F{-,z) 

l/A Jr{l/uj)=r{l/T)=p ^ 



X exp 



1 r^''^ 1 / 



(61) 



where p = AT/uj and the limit AT~^ — )■ oo is taken. The string (surface) tension parameter 
l/2a' is introduced. E§ (Note: Dimension of a' is [Length]^. ) The square-bracket ([■ ■ ■])- 
parts of f l6ip are — ^Area = J Vdetgabd'^x (See (175]) ) where Qab is the induced metric 
on the 4D surface. F{p, z) is defined in fl52l) or f lMl) and shows the field- quantization of 
the bulk scalar (EM) fields. In the above expression, we have followed the path-integral 
formulation of the density matrix (See Feynman's text|35j). The validity of the above 
definition is based on the following points: a) When the weight part (exp [■ ■ ■ ]-part) is 
1, the proposed quantity Scas is equal to E^^^, 0551) . with W = 1 ; h) The leading path 



■^^ a' is a free parameter of the theory. Two typical choices are considered: (a) a' = T ^ (soft surface); 
(b) a' = uj^'^ (rigid surface). 
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is given by rg{z) = l/pg{z), (158|) : c) The proposed definition, (16T|) . clearly shows the 
4D space-coordinates x° or the 4D momentum-coordinates p"" are quantized (quantum- 
statistically, not field-theoretically) with the Euclidean time z and the "area Hamiltonian" 
^ = / V det Qab d^x. Note that F{p,z) or F{l/r,z) appears, in (161 p . as the energy density 
operator in the quantum statistical system of {p°'{z)} or {x°'{z)}. 

In the view of the previous paragraph, the treatment of Sec|6] is an effective action 
approach using the (trial) weight function W{p,z). Note that the integral over (p'^,2;)- 
space, appearing in f lM|) . is the summation over all degrees of freedom of the 5D space(- 
time) points using the "naive" measure d^pdz. An important point is that we have the 
possibility to take another measure for the summation in the case of the higher dimensional 
QFT. We have adopted, in Sec|6l the new measure W{p^^ z)d'^pdz in such a way that the 
Casimir energy does not show physical divergences. We expect the (^zrect evaluation of (16T|) . 
numerically or analytically, leads to the similar result. 



8 Discussion and Conclusion 

The log-divergence in (l54l) is the familiar one in the ordinary QFT. It can be renormalized 
in the following way. 

E^^JAT-^ = -au^ (1 - 4cln(A/w) - 4c'ln(A/T)) = -a(a;,)^ , 

ojr = coi/l- 4cln(A/a;) - 4c'ln(A/T) , (62) 

where Ur is the renormalized warp factor and u is the bare one. No local counterterms are 
necessary. Note that this renormalization relation is exact (not a perturbative result). In 
the familiar case of the 4D renormalizable theories, the coefficients c and c' depend on the 
coupling, but, in the present case, they are pure numbers. It refiects ti^ ^nteract^on 
between (EM) fields and the boundaries. When c and c' are sufficiently small we find the 
renormalization group function for the warp factor u as 

|c| < 1 , |c'| < 1 , Ur = uj{l - c\n{A/uj) - c'ln{A/T)) , 

^(/3-function) = 7— — In — = — c — c' . (63) 
c7(lnA) cu 

We should notice that, in the fiat geometry case, the IR parameter (extra-space size) 
/ is renormalized (see Secjl]). In the present warped case, however, the corresponding 

^"^ In the usual case, the log-terms (divergent terms) are separated and are canceled by the local counter- 
terms. It is difficult, in the present case, to take such a renormalization procedure because the theory is 
free and has no interaction terms (no couplings). The only choice, if we stick to the usual procedure, is 
the renormalization of the wave function and the mass parameter. It does not seem work well. Here we 
take a new approach. We regard the starting boundary parameter w is a bare quantity and the ujr defined 
in ([5^ is a renormalized one. The boundary parameters flow by themselves. No local counter-terms are 
necessary. 

In the list dSS]), all data show |4c| - 10"^ and |4c'| - 10^^ . 
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parameter T is not renormalized, but the warp parameter u is renormalized. Depending 
on the sign of c + c', the 5D bulk curvature u flows as follows. E§ When c + c' > 0, the 
bulk curvature u decreases (increases) as the the measurement energy scale A increases 
(decreases). When c + c' < 0, the flow goes in the opposite way. When c + c' = 0, a; 
does not flow (/3 = 0) and is given by Ur = w(l + cln(a;/T)). 

The final result fl62l) is the new type Casimir energy, —u^. u appears as a boundary 
parameter like T. The familiar one is — in the present context (See ([7]). Note T is the IR 
parameter and is related to / as T = uje~'^K). In ref.|36|. another type T^w^ was predicted 
using a "quasi" Warped model (bulk-boundary theory). 

Through the Casimir energy calculation, in the higher dimension, we find a way to 
quantize the higher dimensional theories within the QFT framework. The quantization 
with respect to the fields (except the gravitational fields Gab{X)) is done in the standard 
way. After this step, the expression has the summation over the 5D space(-time) coordinates 
or momenta J dz Yla dp""- We have proposed that this summation should be replaced by the 
path-integral J Yla z^P'^i^) with the area action (Hamiltonian) A = J ^ det Qabd'^x where 
Qab is the induced metric on the 4D surface. This procedure says the 4D momenta p"" (or 
coordinates x"") are quantum statistical operators and the extra-coordinate z is the inverse 
temperature (Euclidean time). We recall the similar situation occurs in the standard string 
approach. The space-time coordinates obey some uncertainty principle |37j. 

Recently the dark energy ( as well as the dark matter ) in the universe is a hot subject. 
It is well-known that the dominant candidate is the cosmological term. We also know 
the proto-type higher-dimensional theory, that is, the 5D KK theory, has predicted so far 
the divergent cosmological constant [5]. This unpleasant situation has been annoying us 
for a long time. If we apply the present result, the situation drastically improve. The 
cosmological constant A appears as 

1 matter 

S = j d^x^/^{-^{R + \)} + j d'^x^f^{Cmatter} , g = deig^^ , (64) 

where G^v is the Newton's gravitational constant, R is the Riemann scalar curvature. 
We consider here the 3+1 dim Lorentzian space-time (/i, = 0,1,2,3). The constant A 
observationally takes the value. 



^Kbs - ^ \ 2 -<- {l^-''eVY , Ao,, ~ ~ 4 X 10-66(eV^)2 , (65) 



where Rcos ~ 5 x lO^^eV ^ is the cosmological size (Hubble length), m^j is the neutrino 



Both IR and UV boundary interactions (c' and c ) contribute to the scaUng behavior of the system. 
This should be compared with the usual perturbative (w.r.t. the coupling) case where (3 is calculated from 
one region, for example, UV-region. 

This is the case of "asymptotic free" in the usual renormalization of the gauge coupling of 4D YM. 
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mass. On the other hand, we have theoretically so far 

1 . 1 



^ A,, ~ = M/ ~ (lO^Wr . (66) 



This is because the mass scale usually comes from the quantum gravity. (See ref . |40) 
for the derivation using the Coleman- Weinberg mechanism.) We have the famous huge 
discrepancy factor: 

- ■ , Ndl = MpiR,,, ~ 6 X 10^° , (67) 



A 



obs 



where A^^)^ is the Dirac's large number [H]. If we use the present result (162|) . we can obtain 
a natural choice of T, w and A as follows. By identifying T~'^Ecas = — «iAT"^a;^/T^ with 
/ d^Xy/^{l/Gj\f)\ob = Rcxisi^/^N), we obtain the following relation. 

1 C(j^ A 

^DL = ^cos~pr~ — ""^i^T^ ' '^1 • some coefficient . (68) 
Gat 

The warped (AdSs) model predicts the cosmological constant negative, hence we have 
interest only in its absolute value. El We take the following choice for A and u. 



A = Mpi ~ lO^^GeV , Lu ~ ^ = J ^ ^ 10"^eV . (69) 

The choice for A is accepted in that the largest known energy scale is the Planck energy. 
The choice for u comes from the experimental bound for the Newton's gravitational force. 

As shown above, we have the standpoint that the cosmological constant is mainly made 
from the Casimir energy. We do not yet succeed in obtaining the value ai negatively, but 
succeed in obtaining the finiteness of the cosmological constant and its gross absolute 
value. The smallness of the value is naturally explained by the renormalization flow as 
follows. Because we already know the warp parameter u flows (l63l) . the Xobs ~ l/-Rcos 
expression (169|) . Xobs oc w^, says that the smallness of the cosmological constant comes from 
the renormalization flow for the non asymptotic-free case (c + c' < in f l63|) ). El El 

The IR parameter T, the normalization factor A/T in fl54p and the IR cutoff /i = A^ 
are given by 



The relation ~ y/ Mpi /Rcos = \/^/RcosVGn, which appears in some extra dimension model[38l 
139) . is used. The neutrino mass is, at least empirically, located at the geometrical average of two extreme 
ends of the mass scales in the universe. 

This fact strongly suggests de Sitter (dSs) version of the present work could solve this sign problem. 
Details are under way. 

A.M. Polyakov presented, in an early stage, the idea that the cosmological constant may be screened 
by the IR fluctuation of the metric [32]. It was clearly shown, using the 2 dim R^-gravity, such thing really 
occursfiS], 

We claim here the smallness of the cosmological constant is dynamically explained (without fine 
tuning) . 
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where is the nucleon mass. El The FigJT2] strongly suggests that the degree of freedom 
of the universe (space-time) is given by 




~ 10 



,74 




(71) 



In the recent exciting work by Hofava|^ H5]. the vastness of the string theory is stated 
as "the string theory represents a logical completion of quantum field theory, not a single 
theory". He tackles the renormalizability problem of the quantum gravity not from the 
string theory but from a "small" one, that is, a 3+1 dim local field theory with spacially 
higher-derivative interactions. The idea comes from the success of Lifshitz theory |16], a 
spatially-higher-derivative scalar theory, in the condensed matter physics. The present 
approach has some points which can be compared with Hofava's. Basically both are the 
gravity-matter local field theory inspired by the string, brane and membrane theories. 
Hofava's one is basically 3+1 dimensional, while the present one 4+1 or 5 dimensional. 
In both ones, the renormalization flow plays a key role in the renormalizability (UV- 
completion). In the Hofava's, the Lorentz symmetry is abandoned as the starting principle, 
but is regarded as the dynamically emergent one (in IR region of RG flow). At the cost 
of Lorentz symmetry, he introduces spatially higher-derivative terms in order to suppress 
divergences in UV-region. This anisotropy between space and time (non-relativistic aspect) 
does not appear in the present approach because we treat the 4D world isotropically. 
Instead we have anisotropy between the 4D world and the extra space. In the present 
case, the renormalizability is realized by the warped configuration (thickness) and the 
appropriate suppression by the weight function. We need not higher-derivative terms. The 
origin of the suppression is, at present, not established. Also in Hofava's, an uncertain 
procedure "detailed balance condition" is introduced in order to reduce the number of 
independent coupling constants. It looks important to take a new standpoint or view, 
which has been overlooked so far, about some basic things, such as the starting symmetries, 
the quantum treatment of gravitational and matter fields, the regularization method and 
the meaning of the extra axis(es), to solve the divergence problem of gravity. 



9 App. A. Equation of Minimal Surface in AdSs Ge- 
ometry and Classification of Surfaces 



The present new idea is that the regularization surfaces are determined by the principle 
of the minimal surface in the bulk AdSs manifold. We require the ultraviolet and infrared 
regularization surfaces obey the law of the higher dimensional (bulk) geometry. In this 
section, we examine the minimal surface for the warped and fiat cases. We classify all paths 
(solutions). It is useful in drawing minimal surface curves in the text and in confirming 

Note that the present model predicts the nucleoli mass scale (1 Gev) from the 3 data: Planck mass 
Mpi, the neutrino mass m^. and the cosmological size Rcos (or the cosmological constant Xobs ^ Rcos)- 
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the non-crossing |fj of curves. 

The AdSs geometry is described by 



1 



a,b= 1,2,3,4 



1 1 

- < z< - 
u - - T 



(72) 



where 6ab is the 4D Euchdean flat metric (not Minkowski) {6ab) = diag(l, 1, 1, 1). S^-sphere 
in the "3-brane" located on z of the extra coordinate is expressed as 



[X 



1\2 



(X 



2\2 



[X 



3\2 



[X 



4\2 



(73) 



where we allow the radius to change along 2;-axis. r{z) is some function of z which must 
be determined by the minimal area principle in the AdSs geometry. On the 4D surface, in 
the bulk, defined by (17H|) , which we call the regularization surface or boundary surface, the 
line element can be expressed as 



ds^ 



1 



{5 ah + r '^,,^ )dx"-dx'' = gab{x)dx°-dx'' 



uj'^z'^ 



[rr 



dr 
dz 



(74) 



The surface area is given by 



A = j y/det gab d'x = j p^yi + ^c?' 



X 



l/T ^ 



1/. ^'z^ 



1 r^dz 



A 



Ip- 



1/^ u^z^p3]l p^ 



Idz 



p 



(75) 



Under the variation r{z) — )■ r{z) + 6r{z), A changes as 
1 



6A cx 



1^' ijn'-' 



6r 



1 ^/T 




+ 
- 1/ui 


Ji/oj L 



With the condition that the radii at the end-points are fixed: 



5r 



2=l/a; 



5r 



z=l/T 







' \ 6rdz . (76) 



(77) 



we obtain the differential equation of the minimal surface 

, 

z r'^ + 1 
4 p' pp" — 2p"^ 



3 H — r r 



1 1 

- < z< - 
u - - T 



3- 



z p^ 



p'2 _|_ p4 







p 



r 



1 1 

- < z< - 
u - - T 



(7J 



This requirement comes from the renormahzatioii interpretation of the present regularization. See a 
few hnes above (051) • 
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In terms of the coordinate y [uoy = \n{uz) for ojz > 1, y > 0), (ITHj) can be expressed as 

3e'-y + Acorr - ' = , f^^ , r ^ ^ , 0<y<l . (79) 

Let us examine the trajectory r = r(2;) or r = riy) which are the solution of (175]) or ( I79p 
respectively. 

(A) Flat limit 

In the y-expression f l7U|) . we can take the flat limit: w = 0. 

flat limit (w = 0) 



, < y < / . (80) 



1 +f2 

In terms of m = l/r^, the above one can be expressed as 

1 1 

u[y) = , ,^ = > 

r{yY x'^x" 
flat limit(cj = 0) : u = -Qu^ < , < y < / 

From this equation we know an important inequality relation: 

I 

2 



u\y=i — u\y=Q = ~Q u dy < . (82) 
Jo 

The inequality n < in flHTl) implies n(?/) is convex upwards. 

Making use of the above relation, we can classify all solutions as follows. 
(i)u(y = 0) > 

(ia) n(/) > 

In this case u{y) > for < y < I. u{y) is simply increasing (r(y) is simply 
decreasing) , 

FigED] 

(ib) M(/) < 

(iba) w(0) < u{l) , Figi 
(ib/3) u{0) > u{l) , Fig( 
{ii)u{y = 0) < 

u{y) is simply decreasing (r(y) is simply increasing), Sample 7, Figj23] 

Although numerical solutions are displayed in Fig EUll^ the flat limit case is exactly 
solved and was explained in Appendix A of ref.|2n]. We have confirmed the high-precision 
equality between the numerical curves and the analytical ones. 

(B) Warped Case 
In terms oi u= l/r^, eq.( !78|) can be rewritten as 

z u-^ u' + 4n'^ dz 
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Figure 20: Geodesic Curve (IHOl) 
by Runge-Kutta. Type (ia). 
r(0) = 4.472, r(0) = -22.36. 
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Figure 21: Geodesic Curve (IHOl) 
by Runge-Kutta. Type (iba). 
r(0) = 2.236, r(0) = -5.590. 
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Figure 22: Geodesic Curve (IHOl) Figure 23: Geodesic Curve 
by Runge-Kutta. Type (ib/3). ( IHOi) by Runge-Kutta. Type 
r(0) = 1.4142, r(0) = -1.4142. (ii).r(l.O) = 10.0,r(1.0) = 350.0 
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From this equation, we obtain the following inequality relations. 



d . , ,2 ■ _ 1 



i2 



dz 



> 



z u-^ 



n'(6n2 + u" 



1 u'^ 



-dz > 



Z W 



> 



z u-^ 



Note that the second equation implies {u' + Au^)\iit > {u' + 4m^)|i/^j. 
Bl) z-flat limit 

Before the classification of all solutions, we note here, in the case that the sphere radius 
slowly changes (flows) in the following way. 



dr 



dz 



equivalently 



d{r' 



d{z^ 



< 1 



^5) 



the equation flTSj) reduces to the z-flat limit: 



r'2 + 1 



d 

dz 



< z< 



16) 



which is the same as f lHOl) except the variable z@ The classification of the z-fiat solutions 
goes as in the previous case. 
{i)u'{z = l/cu) > 

(ia) u'(z = 1/T) > 

Simply Increasing (r(z) is simply decreasing). 

(ib) M'(^ = 1/T) < 

(iba) u{z = l/uj) < u{z = 1/T) 
(ib/3) u{z = l/tu) > ulz = 1/T) 
{u)u'{z = l/co) < 

Simply Decreasing (r(z) is simply increasing) 

B2) General Case 

Let us consider the general warped case. The first inequality relation of (18^ implies 



u 



i2 



Au is simply increasing for l/a;<z<l/T. 
{S){u'^ + Au')\,/T > 1 

(ia) K' + 4M3)|i/, > 1 

u'^ + 4^3 > 1 for 1/cj < ^ < 1/T. 

(ib) (M'' + 4^3)1,/^ < 1 

a < n'^ + Au^ < (3 for 1/uj < z < 1/T, where a and /3 are the constants defined 
by < a = (n'' + Au^)\i,^ < 1, /? = + Au^)\i,t > 1 
{n){u'' + Au')\,/T<l 



^"^ It is interesting that, in this hmit, the parameter 1/uj appears as the UV-cutofF and 1/T appears as 
the IR-cutofF of z-integraL Compare with the flat geometry case summarized in SecUJ 
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Figure 24: Numerical Solution by Runge-Kutta. ([83]), l/w = 10^^ < ^ < 1.0 = 1/T. 
= 0.9~^,u'(l) = 0.0 . Graph Type (ia). Vertical axis is u = 
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1.2 



1 . 1 



. 9 




1 



u'^ + 4^3 < 1 for l/w < 2 < 1/T. 
Note the relations: -u'^ + 4m'^ = 4(1 + r''^yr^ = zu^lGu"^ + u")/u' = —z{3 + 3r'^ — rr")/r'r'^ 

Taking into account the last inequality relation of (184|) . the above three cases satisfy 
the following relations between u, u' and u" . 

(ia) When n' > is valid, u" > -6u^ , u'^ > 1 - 4u^ 
When u' < is valid, u" < -6n^ , n'^ > 1 - 4n^ 

See FigMM 

(ib) When n' > is valid, u" > -6n^ , /3 - 4n^ > n'^ > a - 4m^ 
When n' < is valid, u" < -6n^ , (3 - 4n^ > n'^ > a - 4u^ 

See FigJ2SlEl 

(ii) When m' > is valid, u" > -Qu^ , m'^ < 1 - 4m^ 
When m' < is valid, u" < — , m'^ < 1 — 4m^ 
See FigJSiEH] 

10 App. B. Casimir Energy of 4D Electromagnetism 

We review the ordinary Casimir energy of the 4D electromagnetism in the way compa- 
rable to the 5D analysis in the text. □ The content is described in a way suggestive to 

See, for example, the text|47|. 



43 



Figure 25: Numerical Solution by Runge-Kutta. ([83]), 1/w = 10^^ < ^ < 1.0 
u{l) = 0.9~^,'u'(l) = 0.0 . Graph Type (ia). Vertical axis is r = 1/a/u 
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Figure 26: Numerical Solution by Runge-Kutta. (ES]), l/u = 10"^ < ^ < 1.0 = 1/T. 
u{l) = 1.0,m'(1) = 0.5 . Graph Type (ib). Vertical axis is m = 1/r^. 
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Figure 27: Numerical Solution by Runge-Kutta. (183|) . l/w 
u{l) = 1.0,m'(1) = 0.5 . Graph Type (ib). Vertical axis is r = 

1.4- 
1.2 
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10^^ < ^ < 1.0 = l/T. 
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Figure 2 
= 
0.8r 
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3: Numerical Solution by Runge-Kutta. ([HS]), l/u = 10"^ < z < 1.0 
,5,m'(1) = 0.5 . Graph Type (ii). Vertical axis is m = 



l/T. 
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Figure 29: Numerical Solution by Runge-Kutta. ([83]), = lO""^ < ^ < 1.0 
= 0.5,u'(l) = 0.5 . Graph Type (ii). Vertical axis is r = l/^/u 
2 



0.2 0.4 0.6 0.8 



l/T. 



corresponding quantities appearing in the text. The Lagrangian is given by 

Am = -\F^.uF^'' = \A,d^d^A' + \{d^A^f + tot. deri. , 

F^u = d^A^ - d^A^ , - oo < x'' < oo , (87) 

where the space-time is flat (Minkowski): (?7^jy) = diag(— 1, 1, 1, 1) and fi,u = 0, 1, 2, 3. This 
theory has the U(l) local gauge symmetry. 

A^-^ A^ + d^A{x) , (88) 

where A(x) is the local gauge parameter. We take Lorentz gauge. 

df.A'' = . (89) 

Then the gauge-fixed Lagrangian is given by 

Cem = CTm - \{^^.A^? - tot. deri. = ]^A^d''A^ , (90) 

where 9^ = d^d'' = -{d/dtf + {d/dxf + {d/dyf + {d/dzf, {x^') = {t,x,y,z). Now we 
take the periodic boundary condition of the electromagnetic field A'^(t,x,y,z) as follows: 
for the X and y coordinates it is periodic with the periodicity 2L, and for the z coordinate 
with 21. The former is for the IR-regularization of the two (x,y)-planes, while the latter is 
the separation-length of the two planes. We consider the case L ^ I. 

Periodic B.C.: x ^ x + 2L , y ^ y + 2L , z ^ z + 21 , 

L>/ . (91) 
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Then A'^ is expanded as 

X, y,z) = (flm^m^nW exp{z^(m^x + ruyy) + ijuz} + c.c. 

= (<.m,nWeXp{za;^,r„^r, -f} + C.C.) , 

t^m^myu = ijrm^, ^TUy, jTl), f= {x,y,z) , (92) 

where Z is the set of all integers, and "c.c." means "complex conjugate". The gauge-fixing 
condition (189|) says 



h lUJrrixmyn ' '^mxiriyn 



where [a^) = {a^, a). We take the polarization vector a perpendicular to the wave-number 
(3D momentum) vector uj. 

Then from the gauge condition (193|) . 

a^rrixmyn — coust. (independent of t) . (95) 
Hence aP is not a dynamical variable. Finally we obtain the action of the system. 

fOO 

S = I dt Lem 1 

oo 



L i*L 1*1 

dx dy dzCsM — tot.div. — div. const. 



— 2 2-^ L ■ ^rrixmyTi O'mxmyn ' Q'mxmynS ) 

mix ,my ,n(i7i 

where a)^^^^„ = {oJnixmynf = {rn^^lY + {myjf + {ujf. a* is the complex conjugate of 
a. In the above action, two out of three components of dmxmyn are independent due to the 
relation ( 1^1]) . 

The system (l96l) is the set of harmonic oscillators with different frequencies Um^myn 
corresponding to the degree of freedom of the 3D continuous space. Let us consider the 
quantum mechanical system of the harmonic oscillator. 



Lnoib, b) = ^b^ - ^Cu^b^ , b=^ , -oo<t<oo . (97) 

In order to to examine the quantum statistical property at the temperature T = we 
use the well-known correspondence: Quantum statistical mechanics in D-dim space versus 
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Euclidean quantum field theory in (D+l)-dimension spacetime, < r < /3. (See the 
bookjlH]. D=0 is the present case. ) 



I I dtLHo{-^,b) 



drLuoi'^M 
dr 



t = —IT 

dh 



b' 



dr 



hij) = b{T + /3) : periodic with periodicity f3 . 
The energy of the harmonic oscillator, Eho, is given byF^ 



(9^ 



E 



HO 



Vb LHo{b',b)\fjexp 



- / dTLHoib',b;u) 



d_ 

'd(3 



W{l3,uj) = / Vb exp 



dTLHo{b',b]u) 



,(101) 



where the path-integral is done over all paths which satisfy the periodic condition 6(0) 
6(/3). Using the periodic property fl98l) . 6(r) is expressed as 



2n7i \ ^ 

sm — — r + c„, 

n>l ^ n>0 



2mi 
COS — — r 



(102) 



The first part is odd for the Z2-parity: r — — r and the second one is even. Then W can 
be clearly defined and is evaluated as 



Y\_ dCn Y\_ 

n>l n>0 



exp 



exp 



n>l ^ 



n>0 



n>l 



2n7r, 



n>0 



(.103) 



We normalize W at a; = (free motion). 



Wi(3,0) 
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In terms of the heat-kernel or the density matrix (See the text by Feynmanf 



< xle-f^^lx' >= p{x, x'; /3) = / X»&(t) exp 

Jb{0)=xMP)=x' 



(104) 



(99) 



(Note that, in text, the periodicity (^5)) requires 6(0) = 6(/3) = a;. ) we can express 



E 



HO 



Tr < x\Hc-^"\x' / da; < xlHc 



X > 



dx < x\e ^"\x > 



(100) 
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Hence, in order to evaluate W, it is sufficient to consider the following quantity. 



- > m — -75 — — ; — = - azm — — — h / az — -— 

2^ c-fy ^J^oo cj^r J -00 e-2--i 

= ^^ + ^^1 cix^^^ . (105) 
We reach the familiar formula of the energy spectrum of the radiation. 

EHo{f3,i:^) = -g^lnW{(3,u) = '^ + -^^ = Eo{u) + E^{f3,u) . (106) 

The ffist term Eo is the zero-point oscillation energy and does not depend on /3, while the 
second one Ei does depend on /3. We realize the summation over the "Kaluza-Klein modes" 
along the r-direction corresponds to the familiar way of the statistical-procedure over the 
canonical ensemble {exp(— i?„/T), En = {■^ + n)Cj \ = 0, 1, 2 ■ ■ ■ }. This simply means the 
equivalence of the statistical system in the equilibrium (at a temperature T = 1//3) and 
the Euclidean field theory with the periodic Euclidean-time (periodicity /3). 

Going back to the energy evaluation of the 4D EM, we obtain, using the results ( I106p . 

EmEM = Ecas + Ei3 , 
Ecas = 2 EQ^Urn^myn) = ^m^niyn , 

B„ = 2 £.(^„.„.„)=2 5: . (107) 

The factor 2, in front of the middle equations above, reflects the degree of freedom of 
the polarization vector a Ecas is the sum of zero-point energy over all frequency 

modes (vacuum energy of the 4D EM). It is Casimir energy. It is that part of the vacuum 
energy which is independent of the coupling and is dependent on the boundaries. Ep gives 
us Stefan-Boltzmann's law. 

E, = 2 y = 2 r '^^^f = (2L)^(2/) r dkV{P, k) 

= ^(2L)^.2/x- , ViP,k) = --^^ , C(4) = - ,(108) 

where u^k) = k and a formula x^/ {e^ — l)dx = s\C{s + l) is used. V{(3, k) is the Planck's 
radiation formula. The behavior of 'P(/3, k) is graphically shown in Figj2j (We see similar 
graphs in the 5D case of the text. The extra axis corresponds to the /3-axis. ) The peak 
curve of the graph is hyperbolic, /3/c=const., in the (/3, /c)-plane (Wien's displacement law). 
We note Ep oc = T^, hence it vanishes for T = 0. And (2L)^(2/) is the volume of the 
region bounded by the two planes. 
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Ecas does not vanish for T = 0. It is, however, formally divergent. We need a proper 
regularization for the summation over the infinite degree of freedom due to the continuity of 
the space-time. It corresponds to the renormalization procedure in the local field theories. 
From the explanation so far, Ecas is given by 



^Cas 



E 



mxiriyn 



9{ 



rux ,my ,nGZ 



rnxniyn % 

A ' 



nix ,my ,nGZ 



•nixniyn % 



(109) 



where we introduce the cni-Oj^ function: g{u/A) = 1 for < cu < A, =0 for cu > A. A is 
the cut-off parameter for the absolute value of the 3D (x,y,z) momentum. We will take the 
limit A — 7- oo at an appropriate stage. We first fix the reference point, L — )■ oo, L ^ I ^ oo, 
from which we "measure" the energy. 



E, 



AO 

Cas 



oo J — oo 



dk~, 



k. 



{Ti/Lf n/iv + ''y + ^^A 

dkxdkydkz /, 2 j_ 1,2 \ t.2 




;iio) 



This quantity diverges quartically. [fj Hence the energy density (the energy per unit area 
of xy-plane) u is given by 



u 



rpA _ 771AO 
^Cas ^Cas 



{2Lf 



oo J —00 



{2ny 



+ + (^7)' 9 {\Jkl + kl + {n 
.nez * V V 



dkz 



k1 + k"^ + kl g y -^^1 k1 + k"^ + kl 



kdk 



k fk 
^9 



vr 

2 -L (n-Y 



dn\ / k"^ 



{n^f 9 [^^lk^ + in^^ 



1 °° POD 

-X(0) + E^W-/ dnXin) 

n=l 



X(n) 







kdkJk^ + in^y g ( 4a/A;2 + (n-^2 



A 



(.111^ 



The last expression of (|110p shows the introduction of the cut-off function g{Ld/A) is equivalent to the 
usual one (fc < A) taken in the text. 
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Using the Euler-MacLaurin formula 



1 f°° 1 1 

-X(0) + X(l) + X(2) + ■■■- dnX{n) = -^B,X'{0) - ^B,X"'iO) + 



B2 = - , ^4 = -— , (112) 
6 30 ^ ^ 

where i?„ is the Bernoulh number, we finally obtain the finite result. 

_ TT^ ^4 _ vr^ 1 

(27)3 IT ~ "720 (203 ' (113) 

which does not depend on A. Especially there remains no logA divergences. This 
point is contrasting with the ordinary renormalization of interacting theories such as 4D 
QED and 4D YM. Hence we need not the renormalization of the wave-function and the 
parameter /. As is shown in the text, the renormalization of the boundary parameter(s) is 
necessary in the 5D case. 

11 App. C. Numerical Confirmation of the Relation be- 
tween Weight Function and Minimal Surface Curve 



In this appendix, we numerically confirm the proposal in SecITl In order to define the 
weight function W{p,z), we presented the requirement (160!) , that is, the valley-bottom 
line of the (p, z)-integral of (I55p should be equal to the minimal surface line For 
the case of the linear suppression (the weight W4, Fig JTTI) . its valley-bottom line is read 
from the contour graph of FigjSOl In FigJSTl the line is numerically reproduced as the 
minimal surface line. For most of other suppression forms, we confirm their valley-bottom 
lines can be reproduced as the minimal surface lines by taking the boundary conditions 
appropriately. 



12 App. D: Numerical Evaluation of Scaling Laws: Ecas 

m, E§i m, and m 

In the text, (regularized) Casimir energy is numerically calculated in three ways: 1) Original 
version (Rectangle- region integral), 2) Restricted-region integral ( Randall-Sundrum type ), 

This means the interaction between (4D) free fields and the boundary is so simple that there is no 
anomalous scaling behavior. 5D free fields, however, turn out to have the anomalous scaling behavior as 
is shown in the text. 

Note that the positive definite expression (|109p is, after the change of the energy origin (jllip . assigned 
to a negative value. The experimentally-observed attractiveness of the Casimir force tells the importance 
of this regularization procedure. 
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Figure 31: Minimal Area Curve p{z), (178|) . p(l.O) = 30.0, p'(l.O) = 10.0. Horizontal axis: 
0.0001 < z < 1.0 ; Vertical axis: 0.0 < ^ < 30. 
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Figure 32: Casimir Energy E^^^ of flll4p for various A. T = 1, w = 10"^. Horizontal axis: 
In A (A = 10^ X (1,2,4,8,16)), Vertical Axis: -\n{\2^TT^ x (1/2)E^;^^|). 
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3) Weighted version. The final expressions show the scaling behaviors about the boundary 
(extra-space) parameters T , u and the 4D momentum cut-off A. The results are crucial 
for the present conclusion. Hence we explain here how the numerical results are obtained. 

First, let us take the un- weighted case with the rectangle integral-region (original form) 
of Casimir energy (l38l) . 

- ^Sa! ^) = ^'^^ ' 

F^{p,z) = j^^j^~kGUz.z)dk . (114) 

where is explicitly given in (133|) . The integral region is graphically shown, in Fig 13 
as the rectangle ([l/w, 1/T] x [/i = AT/a;,A]). The graphs of the integrand of f lll4p . 
(-l/2)p3F(p,2), are shown for (T,w,A) = (1, 10^ 10^)[Fi^IH], (1, 10^ 2xlO^)[Fi^E], (1, 10^ 4x 
10'*)[-Fi(7|TU], in the text. From the behaviors we can expect E^^^{uj,T), ([mD, leadingly 
behaves as A^/T, because the depth of the valley, shown in FiglHlfTOl proportional to A^ 
and their behaviors are monotonous (except near the boundaries z = 1/u and 1/T) along 
the extra axis. It is confirmed by directly evaluating flll4p numerically (the numerical 
integral in |49j). We plot the numerical results in Fig|32]for various A's with fixed (T, w), 
and in FigJ33]for various T's with fixed [u, A). Furthermore we have confirmed sufficient 
w-indepedence for the region T = l,w = 10^ x (1, 2, 4, 8, 16), A = 10"^ x (2,4). Let us fit 
E^^giu^T) from the above numerical results. First we can regard it as the function of 
one massive parameter A, and two massless parameters A/T and K/uo. From the linear 
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Figure 33: Casimir Energy E(j'^^ of (I114p for various T. uj = 10^, A = 2 x 10^. Horizontal 
axis: InT (T = (1,2,4,8,16)), Vertical Axis: \yi{\2^tx'^ x {l/2)E^'~\). 
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dependences in FigJ32] and FigJSSl we may put 

A^ A A 



E^„,(a;,T) = — (ai + a2ln- + a3ln-) • (115) 

From the cu-independence, we take = 0. From the numerical results, the best fit is given 
by (Manipulating Numerical Data in |49| ) 



-0.02500 — (1 -4.685 x 10"^ In-; 



;il6) 



From the precision of the numerical integral, we may safely regard the log term in flll6p 
vanishing (02 = 0). Finally we obtain ( l39i) . 

For the restricted region case, ( 14T|) . we do the numerical integral of the following ex- 
pression. 

- Ecj^fiuj, T) = —ldq dz q'F-{q, z) = -— du dp p'F~{p, u) (,117) 

where /i = AT/u. We plot the results, in FigJMt for various A with fixed {u},T) and in 
FigJ35l for vaious u with fixed (T, A). Furthermore we have confirmed T-independence for 
T = (1, 2, 3) with {u = 10^ A = 2 x 10^) and for T = (1, 2, 3, 4, 8, 16) with {u = 10^, A = 
2 X 10^). From the straight line behaviors and the T-independence, we can safely fit the 
curve as E^^^ = (AVu;)(6i + 62 ln(A/r) + 63 ln(A/u;)), 62 = 0. The best fit is given by 

2^TT'^Ecaf{u,T) = — ((-1.59, -1.56) x 10"^ + (-1.41, -1.97) x 10"^ ln(A/w)) .(118) 

CO 
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Casimir Energy E^j^f of flll7p for various A with fixed 



Figure 34 

Horizontal axis: InA (A = (1,2,4,8,16) x 10^), Vertical Axis: ln{\2^7i'^ 
results are placed on a straight line with the slope 5 for different A's. 
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Figure 35: Casimir Energy E^f/ of ffTTTD for various u with fixed (T = 
Horizontal axis: \nu {u = (1,2,4,8,16) x 10^), Vertical Axis: \n{\2^TT^ 
results are placed on a straight line with the slope -1 for different cu's. 
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Figure 36: Casimir Energy E~Y/ of (fTT9|) for various A with fixed (T = l,u! 
103)(RightDown) and (T = 0.01, w = 102)(LeftUp). Horizontal axis: In A (A 
10^ X (1,2,4,8,16) and 10^ x (1,2,4,8,16)), Vertical Axis: -\n\2^Tr^{^)E-J^'\. 
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The first component of the above coefficients comes from FigJM] data, the second one from 
FigJSSldata. The "width" of the coefficient- values tells us the first-term coefficient has the 
significant digit number 2, while the second-term one has the number 1. In the text, we 
take the average values (HI|. 

Finally we explain the weighted case (152|) taking the elliptic type, Wi, as an example. 

Jp<A [^T^r Ji/uj 

(iV^)-ie-(V2)pV-^-(i/2).^T2 iVi = 1.711/87r2 , (119) 

where the UV cut-off A is introduced to see the scaling behavior. In FigJ5^ we show the 
numerical results of E^^/{uj,T) for different A's with fixed (T, w). The two lines both are 
straight ones with the slope -1. For the T-dependence (with fixed [u, A)) and w-dependence 
(with fixed (T, A)) we show them in Figj37]and in Figj38| respectively. They are straight 
lines with slopes +1 and -4, respectively. From the straight-lines behavior of Fig JHUllHH| we 
can safely fit the curve as E^^^ = (a;^A/T)(ci -|- C2ln(A/a;) + C3ln(A/T)). The best fit is 
given by 

(2ViVi/2) X E^^^(w,T) = ^A (^-1.04-0.11 In ^ + 0.099 In . (120) 

Taking into account the present precision, we take Ci = 1.04 x (2/1.711) = 1.22, C2 = 
0.11 X (2/1.711) = 0.13, C3 = -0.099 x (2/1.711) = -0.12 in the text. As for other types 
of ly's, the best fit scaling behaviors are listed in (153|) of the text. 
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Figure 37: Casimir Energy E^Z' of (ITT9|) for various T with fixed (A = 2 x 10^, u = 10^). 
Horizontal axis: In T (T = (1,1/2,1/4,1/8,1/16)), Vertical Axis: -\n\2^'ir^{^)E-^Y'\- 
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Figure 38: Casimir Energy E(j^^ of flllOp for various uj with fixed (A = 8 x 10^, T = 1). 
Horizontal axis: \n u {u = 10^ x (4, 8, 16, 32, 64, 128)), Vertical Axis: - In \2^TT^{^)EcaY'\- 
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13 App. E. Normalization Constants of Weight Func- 
tions ( 1521) 



In SecJHl we introduce the weight function W to evaluate the Casimir energy. For the 
comparison between the Casimir energy values obtained by different W^s, the normalization 
constants are important. The normalization constants Ni^s fl52|) are defined by the following 
condition: 

dz Wi{p, z) = -z-i^— / dx dw x^Wi{p = cox, z = w/T) 



^ , A>c^>T (^21) 



For the ends of the integral-regions, we practically may take ^ = oo and ^ = except for 
W2 and W^. As for the starting end of x-integral, we have two choices depending on what 
range of the value u is considered ( (A) uj ~ \/ AT, (B) cu ^ VAT (C) u) <^ \/ AT ) in the 
numerical data-taking. They are explicitly given by 

(A) cu ~ -\/AT (geometrically averaged point) 
In this case, we take ^ = 1, ^ = 0, ^ = cx) in flT^ . 

o /•! o roc 

Stt^Ni = — dw e-'"'/2 = 1.557 , Sn^Nib = ^= dx a^^^e'^'/^ = 1.820 
Ve Jq Ve Ji 

Stt^N^ = f dx I dwxh-'''" = 2(^)3 



1 Jt/uj 
00 pi 



Stt^N, = I dx I dwx^'e-^"^"/^ = -{"if 



1 Jt/uj 



POO /•! POO p1 

Sti^N^ = dx dw x^e-^'/^^' = 0.3222 , Stt^N^ = dx dw x^-^^""' = 0.6342 

Ji Jo Ji Jo 

^TT^A^fi = / dx I dw x^e-^'/2"' = 0.5521 , Sn^Nj = / dx x^q-""'/^ ^2^^ ^'^^^^ 



7 

1 JO Jl 



3 

V^Jo 



1 POO pi 

iTT^N^ = A= I dw e"^/"'"' = 0.3800 , Sn'^N^r = j dx j dw a;3e-'(-'+V-^)/2 _ 0.03893 

/oo p1 
dx j dw a;3g-(x/«,)(x+i/«,)/2 _ Q^j34g 

(*oo p1 

8n'N,, = I dx dw a:h-imi-'+y^')' = 0.005006 



1 JO 



/oo pi 
dx J dw a;3e-(V2)(^+i/«')^ = 0.03921 (,122) 

(B) u > Vat 
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In this case, we take ^ = 0, ^ = 0, ^ = oo in (1T2B . 

p1 poo 

STT^iVi = 2 X / dw e-'"'/^ = 1.711 ^ Sn^Nu = dx x^^e'^'/^ = 2 

87r2iV2= I dx I teV^'" = 2(''' " 



JT/w ^ 

87r^N,= I dx C dwx\-^"^"'^-'^''^^^ 



T/uo 3 T 

87r2Ar4 = / dx dw a:=^e-^'/2"'' = - , Svr^iVs = / dx dw x^^e"^/"'' = - 
Jo Jo 5 Jo Jo 3 

/•oo /•! r) /'OO -| 

Stt^Nq = dx dw x^e-^'/^^; _ _ ^ §^2^^^ = dx x^e'^'/^ = - 
Jo io 3 Jo 2 

Stt^Ns = 2 X /" rfw e-^/2"'' = 0.4177 
Jo 

POO P 1 

Sn^N^r = dx dw ^3^-x2(x^+i/«,^)/2 ^ q_^q28 
io io 



/•oo p1 

87r2iV56= / rfx / dw; _ 0.1779 

io Jo 

In'Nss = dx dw a;3g-(i/2)(x2+i/«,2)2 ^ Q_Q^5Q7 
io Jo 

POO p 1 

87r2iV9 = / dx dw a;3g-(i/2)(x+i/«,)2 ^ q Q5320 
Jo Jo 



(123) 



In the text, we take the case (B). 

(C) u < Vat 

In this case, we may take ^ = 0, ^ = 0, ^ = 00 in f ll2ip . This is the same as the case 
(B). In SecJHl we apply the results to the cosmological constant or the dark energy. We 
take there A = lO^^ eV, cu = 10"^ eV and T = lO'^o eV. 
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